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Изменчивость нижней челюсти у лесных мышей, Ѕуѓіуаетиѕ (Мигійае, Койепііа), фауны Украины. 
Лашкова Е. И., Дзеверин И. И., Межжерин С. В. — На примере модельной выборки лесных 
мышей продемонстрирована возможность диагностики 4 видов рода 5улаетих Овпеу фауны 
Украины по мандибулярным признакам. Форма нижней челюсти сходна у разных видов лесных 
мышей, но эти виды хорошо различаются по ее размерам. Точнее других идентифицируется 
$. іаигісиѕ, остальные — менее надежно. 
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Введение 


Нижняя челюсть млекопитающих — информативная и важная для морфологических исследо- 
ваний структура; причем даже у мелких животных зубная, или нижнечелюстная, кость довольно 
прочна и сохраняется в ископаемом состоянии, экскрементах хищников, погадках птиц и т. п. лучше 
других частей скелета. В связи с этим мандибулярные признаки — перспективный объект для коли- 
чественного описания и анализа, как функционального, так и таксономического. В исследованиях, 
посвященых биологии представителей семейства мышиных Мигійае (Койепііа), данные об изменчи- 
вости мандибулярных признаков используют часто и с большим успехом (ЮРетеќег еї а|., 1996; 
Ғеѕііпе, 1972; Кіпрепбеге, Геату, 2001; Ѕсгіуеп, Ваисрам, 1992). 

Цель нашей статьи — продемонстрировать возможности диагностики и количественного 
описания изменчивости по мандибулярным признакам лесных мышей — 5у/улетиѕ Овпеу іп Орпеу 
еі Могобіеу, 1923 — группы грызунов, весьма сложной в таксономическим отношении (Межжерин, 
1997; Миззег еї а|., 1996). Исследование включало в себя общую характеристику изменчивости разме- 
ров и формы нижней челюсти (анализ главных компонент) и проверку возможности определения 
видовой принадлежности (дискриминантный анализ) на модельной выборке. 

Для наглядного представления результатов использованы методы геометрической морфомет- 
рии. Начало данному подходу было положено классической работой Д. Томпсона (Тћотрѕоп, 1942); 
в последние десятилетия он получил широкое распространение и хорошо зарекомендовал себя в 
исследованиях изменчивости формы тела различных организмов (Павлинов, Микешина, 2002 и др.). 
Мандибула мышей весьма удобна для исследований такого рода, поскольку ее форма сравнительно 
проста и может без большой погрешности рассматриваться как двумерная (Оетеќег еї а1., 1996). В 
геометрической морфометрии разработаны методы количественной оценки различий в форме изу- 
чаемых объектов. Мы не используем эти методы в настоящей статье, поскольку основная наша 
задача — характеристика межвидовых различий в мандибулярных признаках — наиболее эффективно 
решается посредством стандартных методов многомерной статистики (Дерябин, 1983). Использо- 
вание нами некоторых приемов геометрической морфометрии позволило сделать полученные нами 
результаты более наглядными. Надеемся, что настоящая статья будет способствовать более широкому 
обсуждению возможностей и ограничений количественных методов в характеристике форм живых 
организмов. 
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Материал и методы 


Материалом для работы послужила коллекция черепов лесных мышей отдела эволюционно- 
генетических основ систематики Института зоологии им. И. И. Шмальгаузена НАН Украины. Были 
получены данные по 89 взрослым особям: 28 экз. малой мыши — 5. иға/епѕіѕ (РаПаѕ, 1811), 15 — 
степной мыши — 5. ағіапиѕ (ВМапЮга, 1881), 18 — лесной мыши — 5. ѕу/уаійсиѕ (Ііппаеџиѕ, 1758) и 
28 — желтогорлой мыши — 5у/уаетиѕ іаигғісиѕ (РаПаѕ, 1811). Схема промеров нижней челюсти была 
заимствована из работы А. Деметра с соавторами (Оетеќег еї а|., 1996). Она включает в себя 13 про- 
меров (СІ, С2, ..., С13) (рис. 1). Фотографические изображения нижних челюстей были получены с 
помощью сканера. 

В работе использованы стандартные статистические методы (Афифи, Эйзен, 1982; Дерябин, 
1983; Клекка, 1989). Общая характеристика изменчивости нижней челюсти получена с помощью ана- 
лиза главных компонент. Собственные векторы и собственные значения вычислены на базе корре- 
ляционной матрицы. Главные компоненты интерпретированы как выражения скрытых факторов, 
определяющих общий характер изменчивости; в анализе использованы три первые главные компо- 
ненты. Рассчитаны и нормализованы факторные нагрузки всех признаков и факторные значения 
всех особей (рис. 2—3). 

Аллометрические зависимости в изменчивости охарактеризованы с помощью коэффициента 
многомерной аллометрии (А), вычисленного для каждого признака у каждого из видов по методике 
П. Жоликёра (ЛоЙсоеиг, 1963, 1984; Дзеверин, Лашкова, 2005). Этот метод также основан на вычисле- 
нии главных компонент, однако для вычисления коэффициентов многомерной аллометрии, исполь- 
зованы ковариационные матрицы лог-трансформированных данных (ЈоЇісоеиг, 1963, 1984) отдельно 
для каждого вида. Вклад 1-й главной компоненты в общую изменчивость оказался заметно меньшим, 
чем при построении главных компонент по корреляционной матрице нетрансформированных данных 
(см. выше), однако во всех четырех случаях превысил 50%. Доверительные границы коэффициента 
были определены посредством бутстрапа особей каждой из выборок (2000 повторностей). Применение 
бутстрапа и интерпретация его результатов основаны на стандартных для биостатистических иссле- 
дований принципах (ЕеІѕепѕѓеіп, 1985). Для проверки предположений о сходстве трендов изменчивости 
использован непараметрический коэффициент корреляции Спирмана (В$): вычислены корреляции 
между значениями коэффициентов аллометрии разных видов по 12 признакам. 

Алгоритмы идентификации видов разработаны на основе линейного дискриминантного 
анализа. Для каждой из моделей определено значение статистики Уилкса (Л). Уровни значимости 
моделей (р) были оценены путем аппроксимации этой величины статистикой Фишера (Е) со сте- 
пенями свободы АН и АР. На основе тех моделей, которые оказались наиболее информативными, 
построены классификационные функции. Методика их использования для определения видовой 
принадлежности лесных мышей подробно описана нами ранее (Лашкова, Дзеверин, 2002). 
Большинство вычислений произведены с помощью компьютерной системы анализа данных 
ЅТАТІЅТІСА (5430, Ірс., 2001, США). Коэффициенты многомерной аллометрии и их доверитель- 
ные границы вычислены с помощью программы РАЅТ (Наттег, Нагрег, 2004). 

Для геометрической морфометрии использованы изображения 18 нижних челюстей (по 5 малой 
и лесной мыши, по 4 — степной и желтогорлой), которые получены с помощью сканера. Для описа- 
ния межвидовых различий в форме черепа использованы 15 меток (рис. 1). Расстановка меток и 
определение их координат произведены с помощью программы ТРЅаіе (Ко, 2001). Вычислены 


Рис. 1. Схема меток и промеров нижней челюсти. 


Ее. 1. Тһе Фаетат оЁ Іапатагкѕ апа теаѕигетепіѕ оЁ тапаіЫе. 
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средние для вида координаты каждой из меток. Эмпирические средневидовые координаты преобра- 
зованы в прокрустовы координаты, по которым построены по методу тонких пластин изображения 
нижних челюстей (рис. 4—7). В качестве эталона использована конфигурация, усредненная для всей 
совокупности данных по четырем видам. Для преобразования эмпирических координат в прокрусто- 
вы и графического представления форм нижней челюсти использована программа РАЅТ (Наттег, 
Нагрег, 2004). 


Изменчивость 


Результаты изучения многомерной изменчивости нижней челюсти лесных 
мышей по малочисленной выборке отчасти сходны с полученными ранее по 
совокупностям экстерьерных и собственно черепных признаков (Лашкова, 
2003 а, 6). На долю первого фактора приходится 87,7% общей дисперсии, на 
долю второго и третьего — соответственно 3,6% и 2,3%. Факторные нагрузки 
приведены в таблице 1. Первая главная компонента может быть интерпретиро- 
вана как статистическая модель фактора общих размеров мандибулы, а осталь- 
ные — как показатели вариаций в форме. Роль этих последних в изменчивости 
нижней челюсти очень низка, и эта изменчивость сводится в первом приближе- 
нии к пропорциональным изменениям. То, что изменчивость в значительной 
степени сводится к различию в размерах, особенно показательно, если учесть, 
что главные компоненты построены на основе корреляционной матрицы, то есть 
по стандартизированным данным (Дерябин, 1983). 

Распределение видов лесных мышей в пространстве мандибулярных главных 
компонент (рис. 2, 3) не отличается приниципиально от такового, полученного 
по экстерьерным и собственно краниометрическим признакам (Лашкова, 2003 а, 
6). Виды лесных мышей заметно отличаются по мандибулярным признакам. 
Лучше всего отграничиваются особи желтогорлой мыши, трех остальных ви- 
дов — значительно хуже. Тем не менее очевидно, что нижнечелюстные призна- 
ки могут быть использованы для видовой диагностики наравне с экстерьерными 
и собственно краниометрическим. 

Таким образом, форма нижней челюсти почти одинакова у разных видов 
лесных мышей, а отличаются эти виды в основном по размерам. Такой тип 
различий характерен для различных групп животных: достаточно часто 
животные разных видов одного рода отличаются по размерам, в то время как 
форма тела животных одного рода остается практически неизменной (Гетеп, 


Таблица 1. Нормализованные нагрузки на первые три фактора, выявленные в изменчивости 
мандибулярных признаков лесных мышей 

Таре 1. М№огтаііхей ѓасёог Іоайіпрѕ оѓ ігѕі @гее ѓас(огѕ ехігасѓей Гот (ће ѕеѓ оѓ уоой тісе тапа шаг сһаг- 
асѓегѕ 


Признак Фактор 1 Фактор 2 Фактор 3 
СІ —0,87 —0,29 0,27 
С2 — 0,94 —0,01 0,15 
Сз —0,82 0,47 0,30 
С4 —0,94 —0,03 0,09 
С5 —0,96 —0,13 0,02 
Сб —0,95 —0,19 —0,07 
С7 —0,97 —0,17 —0,01 
С8 —0,88 0,28 —0,20 
С9 —0,96 0,03 0,00 
С10 —0,95 0,03 —0,13 
СІІ —0,98 0,05 —0,11 
С12 —0,96 0,00 —0,17 


СІЗ 0,98 0,03 0,09 
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Рис. 2. Распределение особей лесных мышей в пространстве 1-го и 2-го факторов: 1 — 5. иға/еиѕіѕ; 
2 — 5. ағіапиѕ; 3— 5. ѕуаіісиѕ; 4— 5. їаигісиѕ. 


Но. 2. ЗсаНегр!оЕ ОҒ Ғасѓог ѕсогеѕ сотршеа Юг ѕресітепѕ оѓ Юпг мооа тісе ѕресіеѕ (РСТ уетзиз РС2): 
1 — 5. иға[епѕіѕ; 2 — 5. ағіапиѕ; 3 – 5. ѕуіуаісиѕ; 4 — 5. іаигісиз. 
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Рис. 3. Распределение особей лесных мышей в пространстве 1-го и 3-го факторов: 1 — 5. иға/еиѕіѕ; 
2 — 5. ағіапиѕ; 3— 5. ѕуаіісиѕ; 4— 5. їаигісиѕ. 


Но. 3. Ѕсаіѓегрої ОҒ Ғасѓог ѕсогеѕ сотршеа Юг ѕресітепѕ оѓ Юпг мооа тісе ѕресіеѕ (РСТ уегѕиѕ РСЗ): 
1 — 5. иғаіепѕі; 2 — 5. ағіапиѕ; 3 — 5. зумайсих; 4 — 5. іаигісиз. 


Егеетап, 1984; Рибо!5, 1988). Представление форм нижней челюсти четырех 
видов лесных мышей как результатов трансформации некоей условной 
усредненной формы (рис. 4—7) наглядно иллюстрирует результаты анализа 
главных компонент. Деформации, необходимые для того, чтобы преобразовать 
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Рис. 4. Обобщенная форма нижней челюсти малой мыши (5. иға/еиѕіѕ). 

Е. 4. СепегаПте4 ѕһаре оѓ тапаібе іп 5. иға/ензіз. 

Рис. 5. Обобщенная форма нижней челюсти степной мыши (5. айапиз). 

Е. 5. СепегаПте4 ѕһаре оѓ тапаіЫе іп 5. апапиз. 

Рис. 6. Обобщенная форма нижней челюсти лесной мыши (5. ѕу/уаїісиѕ). 
Но. 6. Сепега|те4 ѕһаре оѓ тапаіЫе іп 5. зуайсиз. 

Рис. 7. Обобщенная форма нижней челюсти желтогорлой мыши (5. ѓаигісиѕ). 


Но. 7. Сепега|те4 ѕһаре оѓ тапаіЫе іп 5. /аигісиз. 


усредненную форму в ту форму, которая характерна для любого из четырех 
видов, очень незначительны. Решетки, на которые нанесены контуры мандибул, 
почти не искривлены. 

Количественное описание аллометрии в изменчивости мандибулярных 
признаков лесных мышей представлено в таблице 2. Аллометрические диспро- 
порции выявлены для части признаков (для 6 из 13 у 5. игаіепзѕіѕ, одного для 
5. ағіапиѕ, 6 для $. ѕуіуаїйсиѕ и 4 для 5. іаиғісиѕ), остальные признаки можно 
приближенно считать изометричными по отношению к мандибуле в целом. 
Сходные эффекты характерны и для черепа лесных мышей (Лашкова, 2002; 
Дзеверин, Лашкова, 2005). По значениям коэффициентов многомерной алло- 
метрии три вида (5. иға/епѕіѕ, 5. ѕуГуаїісиѕ и 5. іаиғісиѕ) сходны между собой (В$ 
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Таблица 2. Коэффициенты многомерной аллометрии количественных признаков лесных мышей 
четырех видов 


Та е 2. Мшире аЙотейу сое йсіепіѕ юг шапа шаг сһагасќегѕ оѓ юшг жоой тісе ѕресіеѕ 


Вид 
При- 


5. иғаіепѕіѕ 5. ағіапиѕ 5. ѕу/уаїісиѕ 5. іаиғісиѕ 
знак 


А 2,5% 97,5% А 2,5% 97,5% А 2,5% 97,5% А 2,5% 97,5% 
С 3,63 2,44 7,69 100 —156 1,96 5,25 3,72 14,54 1,98 1,01 2,64 
С2 0,44 —0,56 0,94 103 0,54 1,59 1,15 0,41 1,80 1,04 0,66 1,49 
Сз 0,34 —075 2,03 0,94 0,05 1,77 0,13 —2,59 0,85 0,75 0,39 1,27 
Са БИ 0,65 2,10 1,56 1,02 1,84 1,04 0,33 2,74 1,03 0,82 1,29 
С5 114 0,77 146 0,87 0,57 1,40 0,88 2,97 1,66 1,29 0,93 1,58 
с6 БИ 0,52 14 114 0,84 1,64 1,15 0,787 2,01 1,22 0,79 1,58 
С7 114 0,66 1,44 0,85 0,57 1,46 121 0,72 2,02 1,22 0,87 1,57 
С 0,43 —0,28 0,84 0,62 —0,0003 1,08 —0,064 —0,97 0,51 0,97 0,67 1,39 
с9 0,73 0,46 0,92 0,8 0,56 141 0,72 0,62 166 0,56 0,44 0,79 
С10 0,76 0,25 1,03 101 0,75 1,54 0,54 —108 110 0,68 0,41 1,07 
СИ 0,71 —0,12 0,94 107 0,77 136 0,20 —112 0,68 0,81 0,66 0,98 
С12 0,3 —047 1,02 17 0,97 1,61 0,48 —0,66 0,90 0,75 0,005 1,14 
С13 0,73 0,081 0,95 0,98 0,83 122 0,33 —0,61 0,82 0,70 0,52 0,84 


Примечание. 2,5% и 97,5% — нижняя и верхняя границы 95%-ного доверительного интер- 
вала. Достоверные случаи положительной (А > 1) и отрицательной (А < 1) аллометрии выделены 
полужирным шрифтом. 


принимает значения от 0,57 до 0,77; все они статистически достоверны как ми- 
нимум на 5%-ного уровне значимости), а 5. ағіапиѕ заметно отличается от 
любого из них (значения К$ не превышают 0,20, статистически недостоверны ). 
Нет ясности, с чем связано это различие. Это может быть и артефакт, связанный 
с недостаточностью объема выборки, и проявление реально существующих 
особенностей многомерной изменчивости 5. абапи5. 


Видовая диагностика 


Как показали результаты (табл. 3), диагностические возможности признаков 
нижней челюсти лесных мышей не меньше, чем других краниометрических, 
дискриминантный анализ по которым был проведен нами ранее (Лашкова, Дзе- 
верин, 2002; Лашкова и др., 2005). 


Таблица 3. Сравнительная дискриминационная способность разных наборов промеров мандибулы 
взрослых особей четырех видов лесных мышей 


Та е 3. Сошрагайуе @15сгиитапё арііќу оѓ уагіоиѕ сотріехеѕ оѓ тап@ шаг сһагасѓегѕ іп оой тісе 


Доля правильно определенных 
м о Признаки особей А ағ; ағ2 Е М 
50 ЅА $5 УТ | Ср. 
1 13 С1—С1З 96,4 46,7 94,4 100 88,8 0,05 39; 216 9,4 89 
2 9 С2, С4—С6, С8—СИ, СІЗ 96,4 46,7 77,8 100 85,4 0,06 27; 225 14,1 89 
3 7 С2, С4—С6, С8, С10, С12 96,4 40,0 88,9 100 86,5 0,08 21; 2277 15,7 89 
4 6 С2, С5, Сб, С8, С10, С12 96,4 40,0 88,9 96,4 85,4 0,08 18; 226 18,3 89 
5 6 С2, С5, Сб, С8, С10, С12 100 — 88,9 96,4 96,0 0,08 12; 132 292 74 
6 6 С2, С5, Сб, С8, С10, СІЗ 92,9 40,0 94,4 100 86,5 0,07 18; 226 20,5 89 
т 6 С2, С5, Сб, С8, С10, СІЗ 100 — 100 100 100 0,06 12; 132 33,8 74 


Примечание. Уровни значимости во всех случаях: р < 0,0001; М — номер модели; О — 
количество признаков; 50 — 55. иға/епѕіѕ; ЗА — 5. ағіапиѕ; 55 — 5. ѕу/уаїісиѕ; Т — 5. 1аитсиз; Ср. — 
средняя доля. Расшифровка остальных условных обозначений приведена в тексте. 
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Таблица 4. Классификационные функции для определения видовой принадлежности трех видов лесных 
мышей (малой, лесной, желтогорлой) по промерам нижней челюсти (модель 7, табл. 3) 

Таре 4. СазЙсайоп ѓипсііопѕ ѓог @а210$Йс$ оѓ ћгее уоо@ тісе ѕресіеѕ (5. иға/епѕіѕ, 5. зурайсиз, 5. аи- 
ғісиѕ) оп тапд шаг свагасфег$ (шоде 7 гот фаШе 3) 


ИЕ Функция 1 Функция 3 Функция 4 
(55. иғаіепѕіѕ) (5. зу№айси5) (5. ѓаиғісиѕ) 

С2 4,484 6,670 6,839 

С5 0,205 —0,812 —0,035 

Сб —0,799 0,250 —0,509 

С8 2,879 1,816 1,516 

С10 0,567 1,209 1,036 

СІЗ 3,637 3,483 4,420 
Константа — 356,952 —411,530 — 520,568 


Таблица 5. Пример определения видовой принадлежности лесных мышей по промерам нижней челюсти 
(модель 7, табл. 3) 
Та е 5. Ап ехашре оѓ йіарпоѕіісѕ оѓ жоой тісе ѕресіеѕ гот тапііриаг сһагасќегѕ (тойе! 7 гот (аЫе 3) 


Е Реальная видовая Классификационная функция е н ни СОР 
принадлежность 1 3 4 функциям 
1 5. иғаіепѕіѕ 364,38 358,19 345,36 5. иғаіепѕіѕ 
2 5. иғаіепѕіѕ 370,36 359,74 346,45 5. иғаіепѕіѕ 
3 5. іаиғісиѕ 502,40 511,20 516,58 5. іаиғісиѕ 
4 5. ағіапиѕ 378,82 380,92 367,73 5. ѕу/уаїісиѕ 
5 5. іаигісиѕ 485,65 500,96 505,08 5. іаигісиѕ 
6 5. ѕу/уаїісиѕ 408,11 415,01 407,22 5. ѕу/уаїісиѕ 


При идентификации всех четырех видов наилучшей можно считать модель 1 
(табл. 3) с использованием всех 13 признаков. Следует отметить, что, в отличие 
от анализа краниометрических и экстерьерных признаков, когда наилучшим 
образом диагностируются желтогорлая и лесная мыши, по промерам мандибулы 
больше всего правильных определений у той же желтогорлой и — малой мыши. 
Хорошая идентификация желтогорлой мыши, как и во всех других случаях, обес- 
печивается прежде всего за счет больших размеров. 

При анализе только трех видов, за исключением степной мыши (модели 5 и 
7 табл. 3), результаты еще лучше: во второй из них достигается 100% правильное 
определение всех трех видов при использовании всего 6 промеров. Такая 
высокая надежность определения позволяет рекомендовать использование дан- 
ных моделей для установления видовой принадлежности, например, материала 
из погадок, который, как известно, часто бывает представлен нижними челюс- 
тями. В таблице 4 приведены классификационные функции 7-й модели. 

В таблице 5 возможности классификационных функций 7-й модели пока- 
заны на примере идентификации шести произвольных экземпляров лесных мы- 
шей. Ни один из них не был использован ранее, при вычислении классифика- 
ционных функций. Использование классификационных функций позволило 
правильно определить видовую принадлежность 5 экз. Единственный экземпляр 
5. ағіапиѕ был включен в эту группу для проверки возможностей алгоритма. Он 
был определен как представитель весьма сходного с 5. айапиз вида — 5. ѕу/уаїі- 
сиб. Это нельзя считать ошибкой, поскольку возможность идентификации $. ағі- 
апиѕ в 7-й модели не была предусмотрена. 


Мы благодарны И. Я. Павлинову за обсуждение статьи и полезные комментарии, В. М. Титарю 
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